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Abstract. It is reviewed pedagogically how a very successful description of the r\-r\' mass 
matrix can be achieved in the consistently coupled Schwinger-Dyson and Bethe-Salpeter 
approach in spite of the limitations of the ladder approximation. This description is in 
agreement with both phenomenology and lattice results. 



The Schwinger-Dyson and Bethe-Salpeter (SD-BS) approach is the bound-state 
approach that is chirally well-behaved. (For example, see Ref. for recent 

reviews, and references therein for various phenomenological and other appli- 
cations of SD-BS approach.) Therefore, among the bound-state approaches it is 
probably the most suitable one to treat the light pseudoscalar mesons. One solves 
the Schwinger-Dyson (SD) equation for dressed propagators of the light u, d and 
\ s quarks using a phenomenologically successful interaction. These light-quark 

propagators are then employed in consistent solving of Bethe-Salpeter (BS) equa- 
tions for various quark-antiquark (qq) relativistic bound states. Namely, in both 
SD and BS equations we employ the same, ladder approximation, and the same 
interaction due to an effective, dressed gluon exchange. In the chiral limit (and 
close to it), light pseudoscalar (P) meson qq bound states (P = 7t 0,± , K 0,± ,ti) then 
•—. • simultaneously manifest themselves also as (^Mflsi-)Goldstone bosons of dynam- 

ically broken chiral symmetry. This resolves the dichotomy "qq bound state vs. 
Goldstone boson", enabling one to work with the mesons as explicit qq bound 
states (for example, in Refs. [Q,plH,^,g,H,l0[o]) while reproducing (even ana- 
lytically, in the chiral limit) the famous results of the axial anomaly for the light 
pseudoscalar mesons, namely the amplitudes for P — > yy and y* — > P°P + P~ [p"2|] . 
This is unique among the bound state approaches - for example, see Refs. [fl|,|6|,|l2|] 
and references therein. Nevertheless, one keeps the advantage of bound state ap- 
proaches that from the q q substructure one can calculate many important quan- 
tities (such as the pion decay constant f which are just parameters in most of 
other chiral approaches to the light-quark sector. The description [||, |[ ||, |o|] of 
r)-n ' complex is especially noteworthy, as it is very successful in spite of the limi- 
tations of the SD-BS approach in the ladder approximation. 

For the description of r\ and r\ ', the crucial issues are the meson mixing and 
construction of physical meson states. They are formulated in Refs. [^[^,^|] for 

* Talk delivered by D. Klabucar. 
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the SD-BS approach, where solving of appropriate BS equations yields the eigen- 
values of the squared masses, M^ a , M^, and M^ § , of the respective quark- 
antiquark bound states |uu), |dd), |ss) and |us). The last one is simply the kaon, 
and M uS is its mass Mk. Nevertheless, the first three do not correspond to any 
physical pseudoscalar mesons. Thus, M^ Q , M^U, M;r s do not automatically rep- 
resent any physical masses, although the mass matrix (actually, to be precise, the 
non-anomalous part of the mass matrix) is simply Jvl^ A = diag(M^ a , M.^, M 2 g ) 
in the basis qq), (q = u, d, s). However, the flavor SU(3) quark model leads one 
to recouple these states into the familiar octet-singlet basis of the zero-charge sub- 
space of the light unflavored pseudoscalar mesons: 

|7T°) = -L(|uu)-|dd)), (1) 

ln 8 ) = -J=(|uu) + |da)-2|ss)) , (2) 

lno) = -^(|uu) + |da) + |ss)). (3) 

With |uu) and |dd) being practically chiral states as opposed to a significantly 
heavier |ss), Eqs. (|]-^|) do not define the octet and singlet states of the exact SU(3) 
flavor symmetry, but the effective octet and singlet states. However, as pointed out 
by Gilman and Kauffman JT3| ] (following Chanowitz, their Ref . [8]), in spite of this 
flavor symmetry breaking by the s quark, these equations still implicitly assume 
nonet symmetry in the sense that the same states |qq) (q = u, d, s) appear in both 
the octet member n, 8 (||) and the singlet no @. Nevertheless, in order to avoid the 
Ua(1) problem, this symmetry must ultimately be broken by gluon anomaly at 
least at the level of the masses of pseudoscalar mesons. 

In the basis (0-^)), the non-anomalous part of the (squared-)mass matrix of 7t° 
and etas is 

/M 2 n \ 
^naH Ml 8 M| . (4) 
V M§ 8 M 2 J 

The r|8 "mass" Mss = Mr| 8 can be related to the kaon mass through the Gell- 
Mann-Okubo (GMO) relation, although the kaon does not appear in this scheme 
as it obviously cannot mix with n° and etas, since it is strangely flavored. Equa- 
tion (Q) shows that also the isospin I = 1 state n° decouples from any mixing with 
the 1 = states ng and r\o, thanks to our working in the isospin limit throughout. 
Therefore, we are concerned only with the diagonalization of the 2x2 submatrix 
in the subspace of etas in order to find their physical masses and correspond- 
ing qq content. In the isospin limit, obviously M uQ = M d( j, which we then can 
strictly identify with our model n° mass M n . Since in this model we can also cal- 
culate Mj g = (ss|K/l^j A |ss), this gives us our calculated entries in the mass matrix: 

Mfg ^ (n 8 |Ivl 2 NA |n 8 ) = M^ 8 = l[U 2 tt + ~M 2 n ) , (5) 
Mlo ee (rislA^ho) = = ^(M 2 n - M 2 sS ) , (6) 
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M 2 00 ee (no|J&2 A |rio) = + M i) • (7) 

The last one, M 0/ is the non-anomalous part of the r\ "mass" M Tl0 . Namely, all the 
model masses M q q/ (q, q' = u, d, s) and corresponding qq' bound state ampli- 
tudes are obtained in the ladder approximation, and thus (irrespective of what 
the concrete model could be) with an interaction kernel which cannot possibly 
capture the effects of gluon anomaly. Fortunately, the large N c expansion indi- 
cates that the leading approximation in that expansion describes the bulk of main 
features of QCD. The gluon anomaly is suppressed as 1/N C and is viewed as a 
perturbation in the large N c expansion. It is coupled only to the singlet combi- 
nation r|o ([3]); only the r|o mass receives, from the gluon anomaly, a contribution 
which, unlike quasi-Goldstone masses M q q >'s comprising 1vl^ a , does not vanish 
in the chiral limit. As discussed in detail in Sec. V of Ref . in the present bound- 
state context it is thus best to adopt the standard way (see, e.g., Refs. [p^,|l5|]) to 
parameterize the anomaly effect. We thus break the Ua ( 1 ) symmetry, and avoid the 
U a ( 1 ) problem, by shifting the r\ (squared) mass by an amount denoted by 3(3 (in 
the notation of Refs. ||^). The complete mass matrix is then Ivl 2 = M^a + ^A' 
where 

/0 \ 

Ivl A = , (8) 
\0 3(3/ 

and where the value of the anomalous no mass shift 3(3 is related to the topo- 
logical susceptibility of the vacuum, but in the present approach must be treated 
as a parameter to be determined outside of our bound-state model, i.e., fixed by 
phenomenology or taken from the lattice calculations l |l6| ] . 

We could now go straight to the nonstrange-strange (NS-S) basis, but before 
doing this, it may be instructive to rewrite for a moment the matrix (^) in the 
flavor, qq) basis, where 

Ivl A = (3 fl 1 1 j , (9) 

since for some readers it may be the best place to introduce the effect of flavor 
symmetry breaking. Namely, Eq. (g) tells us that due to the gluon anomaly, there 
are transitions |qq) — » q'q'); q, q' = u, d, s. However, the amplitudes for the 
transition from, and into, light uu and dd pairs need not be the same as those 
for the significantly more massive ss. The modification of the anomalous mass 
matrix (||) which allows for possible effects of the breaking of the SU(3) flavor 
symmetry is then 

1 1 X 

iC& = p I 1 1 x I . (io) 



XXX 



2 



There are arguments [|[||], supported by phenomenology, that the transition sup- 
pression is estimated well by the nonstrange-to-strange ratio of respective con- 
stituent masses, X w M. u /M. s , or, as commented below, of respective decay con- 
stants, X f» fn/fss- 
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After adding the anomalous contribution @ to Eq. @, pion still remains 
decoupled and we obviously still can restrict ourselves to 2 x 2 submatrix in the 
subspace of etas. However, when dealing with quark degrees of freedom when 
the symmetry between the nonstrange (NS) and strange (S) sectors is broken as 
described above, the most suitable basis for that subspace is the so-called NS-S 
basis: 

hlN5> = ^duu) + ldd» = -)= |ri 8 ) + y|ho) , (11) 
fl 1 

his) =|ss)=-y-|Tl8) + -^|Tlo) • (12) 

The T|-r|' mass matrix in this basis is 

^2 /M^ ns M* \ /M£ fl + 2P v^X \ >K0 \ 

(13) 

where the indicated diagonalization is given by the NS-S mixing relations)] 

hi) = cos 4>hnv S ) - sin c|j]tis> , hi') = sin 4>\r\Ns) + cos <\>\r\ s ) , (14) 

rotating njvs , r\s to the mass eigenstates r\,r\'. Now the NS-S mass matrix (|3|) tells 
us that due to the gluon anomaly, there are transitions \t]us) his)- However, the 
amplitude for the transition from, and into, t\ns> need not be the same as those 
for the more massive r|s. The role of the flavor-symmetry-breaking factor X is to 
allow for that possibility. As remarked little earlier, there are arguments [||^|]/ 
supported by phenomenology, that the transition suppression is estimated well 
by the nonstrange-to-strange ratio of respective quark constituent masses, A1 u 
and M. s . Due to the Goldberger-Treiman relation, this ratio is essentially equal 
[||||||] to the ratio of njvs and r\s pseudoscalar decay constants f t, ns = f n and f Ms = 
f s§ which are calculable in the SD-BS approach. Same as M u and M s , they were 
found in our earlier papers, especially [^,^,|^]. In other words, we can estimate the 
flavor-symmetry-breaking suppression factor as X w M. u /M. s , or equivalently, as 
X w fn/fss- Our model results M-u./M s — 0.622 and i n /i sS — 0.689 are in both 
cases reasonably close to X exp w 0.78 extracted phenomenologically ||,||] from 
the empirical mass matrix ih^ p featuring experimental pion and kaon masses, 
or, after diagonalization, n and r\' masses - see Eq. (7) in Ref. [^]. (In our model 
calculations below, we use X = f^/fsg = 0.689 to show the robustness of our 
approach to slight variations. Namely, our earlier works lf|[^[l^] mostly presented 
results based on slightly different values of X obtained with the help of ratios of 
yy amplitudes, which is yet another, but again related way of obtaining X.) 

In our present notation, capital M Q 's denote the calculated, model pseu- 
doscalar masses, whereas lowercase ra Q 's denote the corresponding empirical 

1 The effective-singlet-octet mixing angle 9, defined by analogous relations where n N s — > 
ri8,T|s — > rio, 4> — > 6/ is related to the NSS mixing angle tj> as 9 = t|) — arctan \fl — 
4> — 54.74°. The relation between our approach and the two-mixing-angle scheme is 
clarified in the Appendix of Ref. [[|- 
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masses. The empirical mass matrix frl 2 xp can be obtained from the calculated, 
model one in Eq. ( |l3| ) by i) obvious substitutions M uQ = M n — > m n and M s§ — > 
m S s, and zz) by noting that m S s, the "empirical" mass of the unphysical ss pseu- 
doscalar bound state, is given in terms of masses of physical particles as m 2 5 = 
2ra^ — m 2 due to the Gell-Mann-Oakes-Renner (GMOR) relation. Since M utl , ob- 
tained by solving the BS equation, is identical to our model pion mass M w it was 
fitted to the empirical pion mass ra,, e.g., in Ref. Similarly, M u§ = Mk is fitted 
to the empirical kaon mass ttik ■ Therefore we also have M 2 S ss — m 2 , since 
our model has good chiral behavior and also satisfies the GMOR relation, thanks 
to which we have = 2M 2 - — M 2 in a very good approximation. We thus 
see that in our model mass matrix, the parts stemming from its non-anomalous 
part (||) are already close to the corresponding parts in ih^ p . We can thus 

expect a good overall description of the masses in n and n ' complex. We now 
proceed to verify this expectation. 

The anomalous entry (3 is fixed phenomenologically to be |3 exp ~ 0.28 GeV 2 , 
along with X exp ~ 0.78, by requiring that trace and determinant of frt 2 xp have 
their experimental values. But, this can be done also with our calculated, model 
mass matrix )vl 2 . Requiring that the empirical value of the trace m. 2 + m. 2 , ~ 1 .22 
GeV 2 be fixed by Eq. (||), yields 

p = 2Tx 2[(m ^ + m ^' )exp_(M ' Q + M ' s)] (15) 

where X = i n /i sS = 0.689 and lVl ua = 0.1373 GeV and M s§ = 0.7007 GeV 
are now our model-calculated Q quantities, giving us |3 = 0.286 GeV 2 . Since 
( M 2 a + M 2 § ) = 2m 2 c holds to a very good approximation, our approach satisfies 
well the first equality (from the matrix trace) in 

2p+pX 2 =m 2 +m 2 ,-2m 2 , = ^x, (16) 



where the second equality is the Witten-Veneziano (WV) formula [[180, with x 
being the topological susceptibility of the pure Yang-Mills gauge theory. Our 
model values of X and |3 (f^ is fitted to its experimental value) thus imply x = 
(178MeV) 4 , in excellent agreement with the lattice result x = (175 ± 5MeV) 4 of 
Alles etal. @]. 

The mixing angle is then determined to be <p — 43.2° (or equivalently, = 
— 1 1 .5°), for example through the relation 

tan2 * = M 2 — M 2 ' (17) 



where 



T|NS 

and 



M 2 = M 2 a + 2|3 = M 2 +2|3 = 0.592 GeV 2 = (769 MeV) 2 (18) 



M 2 s = M 2 , + |3X 2 = 0.627 GeV 2 = (792 MeV) 2 (19) 
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are our calculated tins and r\s masses. They have reasonable values, in a good 
agreement with, e.g., t|ns and r|s masses calculated in the dynamical SU(3) linear 



a model [17| 



The diagonalization of the NS-S mass matrix gives us the n and r\ ' masses: 

M 2 = cos 2 4> M 2 ns - \/2(3Xsin24) + sin 2 $ M 2 s (20) 
M 2 , =sin 2 4)M 2 NS +%/2|3Xsin24) + cos 2 4)M 2 s . (21) 

Plugging in the above predictions for |3, X, M^,; and M^, our model r\ and r\' 
masses then turn out to be M n = 575 MeV and M^* = 943 MeV. This is in good 
agreement with the respective empirical values of 547 MeV and 958 MeV. 

However, the above is not all that can be said about agreement with experi- 
ment and other approaches. The second thing we may point out is the reasonable 
agreement we find if we insert our values of |3, X and M q q i 's into our model mass 
matrix and compare it with the n-n ' mass matrix obtained on lattice by UKQCD 
collaboration [fjjjj. 

Third, Ref7[H] clearly shows that our approach and results are not in conflict, 
but in fact agree very well with results in the two-mixing-angle scheme (reviewed 
and discussed in, e.g, Ref. [|2(]]). Actually, our results can also be given in 
the two-mixing-angle scheme. 

Fourth, what we found from the mass matrix is consistent with what we 
found in the same SD-BS approach through another route, i.e. from r\,r\' — > yy 
processes [| 1^0] . 

The above shows that the consistently coupled SD-BS approach provides a 
surprisingly satisfactory description of the r\-r\ ' complex, especially if one recalls 
that (3, parameterizing the anomalous r\o mass shift, was the only new parameter. 
Namely, all other model parameters were fixed already by Ref. [£2|] providing the 
model we used in Refs. [|[|^^||^[io]]. Nevertheless, we would like to point 
out that even this one parameter, |3, can be fixed beforehand. Instead of being a 
parameter, |3 can be obtained through WV formula ( |l6| ) from the lattice results on 
the topological susceptibility x- The central value of widely accepted x = ( 1 75 ± 
5 MeV) 4 01 would lead to (3 less than 7% below our model value, which is within 
upper error bar anyway. Thus, we can eliminate (3 as a free parameter and still 
achieve almost as satisfactory description as the one given above. 

Acknowledgment: D. Klabucar thanks the organizers, M. Rosina, B. Golli and S. 
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